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In order to study the mechanism of the mass enhancement in heavy fermion compounds in 
the presence of magnetic field, we study the periodic Anderson model using the fluctuation 
exchange approximation. The resulting value of the mass enhancement factor z _1 can become 
up to 10, which is significantly larger than that in the single-band Hubbard model. We show 
that the difference between the magnitude of the mass enhancement factor of up spin (minority 
spin) electrons z-j -1 and that of down spin (majority spin) electrons zj" 1 increases by the applied 
magnetic field S |[ which is consistent with de Haas- van Alphen measurements for CeCoIns, 
CeRu2Si2 and CePd2Si2. We predict that z^ 1 > z^ 1 in many Ce compounds, whereas z^ 1 < 
z^ 1 in Yb compounds. 
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Heavy fermion compounds such as CeCoIns 1 and 
CeRhIn5 2 have attracted considerable attention because 
of their non- Fermi-liquid like electronic properties, where 
these materials have been considered to locate close 
to the antiferromagnetic (AF) quantum critical point 
(QCP). The anomalous transport properties specific to 
the AF QCP are reported by recent experiments of 
CeCoIn.5. In the normal state of CeCoIns, it is clarified 
that the temperature dependence of the resistivity 3 p and 
Hall coefficient 3 Rh follows p oc T and Rh ocl/T below 
30K. These behaviors are significantly different from the 
Fermi-liquid behaviors, i.e., p oc T 2 and Rjj oc T . From 
a theoretical view point, we have previously derived that 
p oc T and Rh oc 1/T near the AF QCP by considering 
the current vertex corrections. 4 ' 5 

It has been also reported that these anomalous elec- 
tronic properties near the AF QCP are sensitive to a 
magnetic field. For example, p, Rh and v in CeCoIns 
show a nonlinear dependence of the applied magnetic 
field. 3 ' 6 ' 7 Moreover, it is obtained by de Haas- van Alphen 
(dHvA) oscillation measurements that the ratio of the 

spin-dependent effective mass -4 reaches 3 at 15T. 8 ' 9 
Therefore, to understand these interesting phenomena, 
theoretical works about heavy fermion compounds in the 
presence of the magnetic field are needed. In particular, 
it is important to clarify a field dependence of the mass 
enhancement. 

In CeCoIn5, the Fermi surface geometry and renor- 
malized quasiparticle mass have been obtained by dHvA 
measurement. 8 In this compound 4/-electrons mainly 
contribute to the density of state (DOS) at the Fermi 
level. It is known that the applied magnetic field lifts 
the spin degeneracy, and splits the Fermi surface into 
majority-spin and minority-spin surfaces. It is obtained 
that effective masses m* for branch /3 2 of 14th band in 
CeCoIns at 15T are m* = 46m e for the majority spin and 
m* = 154m e for the minority spin, respectively, where 
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m e denotes the mass of a bare electron. 9 The band mass 
for the /3 2 -band is nit, — 1.7m e according to the band cal- 
culation. 8 Then, the obtained mass enhancement factor 
z~ x = m* /rrib is 27 for majority spin and 90 for minority 
spin, respectively. Such a spin-dependent mass enhance- 
ment is widely observed in heavy fermion compounds 
such as CeRu 2 Si 2 , 10 CePd 2 Si 2 n and Celn 3 . 12 

Up to now, electronic properties of heavy fermion sys- 
tems under magnetic field have not been studied enough 
from a microscopic theoretical approach. In the previous 
studies for the impurity Anderson model with particle- 
hole symmetry, it was shown that mass enhancement 
factor (z^ 1 — z^ 1 ) decreases with magnetic field. 13 In- 
troducing the electron-magnon interaction, Edwards and 
Green studied the periodic Anderson model (PAM) in 
strong magnetic field regime, where number of up-spin 
/-electron (nl) is zero. 14 However, dHvA measurements 
are usually performed under a moderate magnetic field 
(~ 10T), where n^/n{ < 1 is satisfied. On the other 
hand, Korbel et al. studied the Hubbard model using 
the Gutzwiller approximation and obtained the relation- 



ship 



l — ni 



for U = 00. 15 Based on this relation- 



ship, — 3r ~ 1.2 is satisfied for — ~ 0.9, which is re- 
alized under 10T according to the magnetization mea- 
surement. 16 However, a dHvA 9 measurement in CeCoIns 
reports ~ 3 under 10T. Thus, it is highly desir- 
able to study the origin of the mass enhancement fac- 
tor by taking the strong spin fluctuation near the AF 
QCP, which is ignored in the Gutzwiller approximation. 
Although Held et al. 17 and Sakurazawa et al. 18 studied 
the Hubbard model under magnetic field using the dy- 
namical mean field (DMFT) method and the fluctuation 
exchange (FLEX) approximation, respectively, they did 
not report spin-dependent mass enhancement. 

The purpose of this paper is to study spin-dependent 
mass enhancements under magnetic field in the PAM. We 
study the two-dimensional (2D) and three-dimensional 
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(3D) PAM based on the FLEX approximation. In the 
PAM, the obtained mass enhancement factor is much 
larger than that in the Hubbard model, because of the 
hybridization between /-electrons and conduction elec- 
trons. In both 2D and 3D systems, the mass enhance- 
ments of up and down spin quasi-particles differ con- 
siderably under magnetic field, which is consistent with 
dHvA measurements in CeCoIn 5 . 9 We also analyze the 
mass enhancement factors using the DMFT in the strong 
correlation regime, and derived a general relationship 
— ~ Pt 1"s , where z„ and £> CT (0) are renormalization fac- 
tor and density of states at Fermi energy with a =\ (j). 
We start with the PAM, 
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[(ef + Ba)flj ka + (e% + Ba)c{ a c ka 



V{!{ a C k(y +c\j ka ) +— ^ /fc T /fc+qT/fc'|/fe'-W) 
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e k = —2ti(cos(k x ) + cos(ky)) — 4t 2 cos(k x ) cos(k y ) 

-2i 3 (cos(2fe x ) + cos(2fe y )) - 2t z cos(fe z ), (2) 

where f ka (f kcr ) is an annihilation (creation) operator 
for an /-electron, and c ka (ct ) is that for a conduc- 
tion electron, respectively, t' and e k are the energy of an 
/-electron and the dispersion of a conduction electron, 
respectively. Ba represents the Zeeman energy, where B 
is magnetic field and a = 1(— 1) corresponds to the T-(-L-) 
spin state. Here, we take ^-factor g and Bohr magneton 
Hb as unity (gfis = 1)- We consider a stacked square lat- 
tice with the Coulomb repulsion U for /-electrons, where 
the hybridization is denoted by V and the intralayer hop- 
ping t\, t2, ts and the interlayer hopping t z for conduc- 
tion electrons. For t z = 0, the above model is reduced to 
the 2D PAM. 

Hereafter, we take t\ = 1 as a unit of energy. We 
apply the FLEX approximation 19-22 where the Green's 
function, the self-energy and the susceptibility are ob- 
tained self-consistently. The FLEX approximation be- 
longs to the "conserving approximations" formulated by 
Baym and Kadanoff. 23,24 From Dyson equation, we ob- 
tain Green's function G^(k) for the /-electron and G c a (k) 
for the conduction electron, 



Gl(k) 
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respectively, where Y> a (k) is the self-energy of /- 
electrons. Here and hereafter, k = (ie n ,k) where e„ = 
(2n + \)nT is the Matsubara frequency. The spin suscep- 
tibilities Xa,a'{q) and the self-energy Yj a (k) in the FLEX 
approximation under magnetic field are formulated in 
Ref. 18 



Mass enhance factor z. is 



_! _ <9Re£ CT (w,fc) 
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where the self-energy T, a (uj,k) represented in the real- 
frequency is obtained by the analytic continuation of 
S CT (ie„,A;) using the Pade approximation. It should be 
noted that the relationship V 2 /(fi - e% F - Ba) 2 <; 1 
is satisfied in actual heavy fermion compounds. In the 
following, we take N = N x X N y X N z = 64 x 64 x 64 
fc-point meshes and the Matsubara frequencies e n from 
-(2N C - l)vrT to (2N C - i)t:T with N c = 512, respec- 
tively. 
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Fig. 1. (color online) Fermi surface in the presence of the mag- 
netic field B = (left panel), B = 0.04 (center) and B = 0.1 
(right) for T = 0.01, (7 = 2 and V = 3. The Fermi surfaces for 
up-spin and down-spin are depicted by the blue (dark) and the 
yellow (light) sheets, respectively. 
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Fig. 2. Mass enhancement factor z„ 1 for an /-electron against B 
for 3D system t z = 0.8 (solid line) and for 2D system t z = 
(dotted line) in the PAM with U = 2, V = 3 and T = 0.01. 
z,7 1 of up (down) spin in the 3D Hubbard model with the same 
Stoncr factor (= 0.96) for B = are depicted by solid (dotted) 
lines. They are almost same lines. 



In the following, we choose the value of the filling 
of /-electrons and conduction electrons as n' = CL 



(n 



Tii), and n c = 0.1 (n c = nS + nf), respec- 



tively, by adjusting /i and e* . The total electron density 
n c . Here, the up(down) spin corresponds 



is n 



= nJ 



to minority(majority) spin. In the actual Ce compounds 
< 1 is satisfied. We put the hopping integrals as 
h = 1, t 2 = -1/6, t 3 = 1/5 and t z = 0.8{t z = 0) for 
three-dimensional (two-dimensional) case. Hereafter in 
the PAM, we fix parameters U = 2 and V — 3. Strong 
AF fluctuation with Q — (it, ir, ir) (3D) and Q = (ir, it) 
(2D) is realized for these parameters. In Fig. 1, we show 
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Fig. 3. DOS of /-electron p a against w in the 3D PAM for B = 
(solid line) and for B = 0.1 (dotted line). The lowest solid line 
denotes DOS for B = on the 3D Hubbard model. 



the calculated Fermi surface in the PAM under the mag- 
netic field B = 0, B = 0.04 and B = 0.1 for T = 0.01, 
where the yellow (light) and blue (dark) sheets denote 
the Fermi surfaces for down and up spin electrons, re- 
spectively. This Fermi surface resembles that of the 14th 
band in CeCoIn 5 . 25 

The mass enhancement factor z^ 1 for /-electrons, 
which is the averaged value of (zfco-) -1 over the Fermi 
surface, is shown in Fig. 2. z^ 1 of up- and down-spin 
electrons at T — 0.01 in the 3D and 2D PAMs are shown 
by full lines and dotted lines, respectively. In the 2D 
and 3D PAMs, we see that the difference between 
and z^ 1 increases with B. It is also shown that z^ 1 of 
the minority-spin {]) electrons is larger than that of the 
majority-spin (J,) electrons. These results in the PAM 
are consistent with the experimental data of dHvA 9 in 
CeCoIns. In Fig. 2, we also plot the mass enhancement 
factors for up (solid line) and down (dotted line) spin 
electrons, respectively, in the 3D Hubbard model. In this 
model, spin dependence of mass enhancement factors is 
not visible (z^ 1 ~ ^j/ 1 )- I n order to realize the same value 
of the Storner factor (= 0.96) in the 3D Hubbard model, 
we choose the filling n = 0.9 and U — 4.6. Thus, the dis- 
tance from the AF QCP is considered to be the same in 
both models. The resulting values of z~ l in the Hubbard 
model are less than 3 and they are much smaller than 
those in the PAM (~ 10). 

In the PAM, the band width without renormalization 
W is estimated by W ~ ~ 2, where D ~ 10 is the 
band width of the conduction band. Then, U/W <~ 1 and 
T/W ~ 0.005 are satisfied in the PAM. We also calculate 
the mass enhancement in the Hubbard model with the 
same U/W ~ 1 and T/W ~ 0.005 as in the PAM. We 
obtain z^ 1 ~ 4.5 for W ~ 8, U = 8 and T = 0.04 in the 
Hubbard model, which is smaller than z^ 1 ~ 10 in the 
PAM. 

Here, the magnitude of the magnetic field B = 0.1 in 
the PAM corresponds to the situation where the mag- 
netic field H ~ 35T is applied to electrons of the 14th 
band of CeCoIns. The reason is as follows. The W in 
CeCoIn 5 estimated from the LDA is - 1000K. The Zee- 
man energy for Ce 3+ is given by (6/7)M[jlbH, where 
6/7 is the g- value and M — 5/2 is spin of Kramers 



doublet. 18 Then, we can estimate the magnetic field 
H = x ±220R _ 35T for B = .1. 

To understand the detailed behavior of the mass en- 
hancement mechanism, we show the DOS p a {ui) of /- 
electrons, which is defined by p a (w) = X)fcA>-fc(w) = 
X^fe -^^- m G a {w + i5, k). In Fig. 3, we see that DOS in the 
PAM is larger than that in the Hubbard model with the 
same Storner factor (= 0.96). The resulting DOS of down 
(majority) spin electrons is larger than that of up (minor- 
ity) spin electrons around u> = (corresponding to Fermi 
energy) under the magnetic field B = 0.1. These results 
originate from the fact that the quasi-particle dumping 
(|ImS|) at the Fermi energy for up spin is larger than 
that for down spin. 

In order to understand the numerical results given by 
the FLEX approximation, we analyze the relationship 
between the self-energy and the resulting DOS based on 
the infinite-order perturbation theory. For this purpose, 
we employ the DMFT scheme where momentum depen- 
dence of the self-energy is neglected. Both in the PAM 
and the Hubbard model, up spin electrons interact only 
with down spin ones. Then, the resulting self-energy can 
be expressed as, 

oo 

S ^) = EE E U n+1 a{T n+1 ,n+l,e v {e^,{e' t }) 

«=o r„ { £l },{<} 

xG CT (iei) • • • G a (ie n )G- a (ie[) ■ ■ ■ G- a (ie' n+1 ), 

where T n denotes a Feynman diagram of the nth or- 
der term, a(T n , n, €j, {ti}, {e-}) is a coefficient of the 
nth order term and {ei}({e^}) is defined as series of 
ei with <t(— a) spin, respectively For example, the co- 
efficient of the second order is a(T2, 2, €j, {e^}) — 
-S(e' 2 - ej +ei - ei). 

Using the spectral representation of the Green's func- 
tion, 

G^) = [ du^±, (7) 
the self-energy is given by, 

oo ( . 

So-^O^EE E doj\ ■ ■ ■ duj n duj[ ■ ■ ■ dco' n+1 

n=0 r n { ei }, {e'jJ 

(7" +1 a(r„ +1 ,n+l,e 3 ,{ ei },{^}) 

x PajuJi) ■ ■ ■ P(y {uj n ) p- a {u[)---p- a {u' n+1 ) 

(ie ni - wi) • • ■ (ie„ n - w n ) (ie n / - wj) • • ■ («e„^ +1 - u' n+1 ) ' 

Here, we assume that the contribution around the Fermi 
energy (cjj = 0) is dominant in the integral for Wj. Then, 
the self-energy can be approximately expressed as, 

oo „ 

S t(* £ j)~EE E / dUJl ' ' ' dbJ ndu'l ■ ■ ■ dw'n+1 
«=0 r n {e t },{<} 

U n+1 a(T n+1 , n + 1, e,, {ej, {<}) JJ(9(n - | Wi |)0(n - M|) 

i 

x Pa(0) n P-.(0)" +1 

(ie ni - wi) • • ■ (ie nn - LO n ) (ie ni - uj[) ■ ■ ■ {ie n ' n+1 - u' n+1 ) ' 
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where il ~ W is a small cut-off energy. According to Eq. 
(9), we obtain 

E T (w)p T (0)~S i Mp i (0). (10) 
Therefore, the following relationship is obtained: 

£^ PiW l ImS T(0)l fnl 
z- 1 ~ p T (0) ~ |ImE,(0)|- 1 ' 

Since all diagrams are taken into account in this anal- 
ysis, the above relationship (11) is expected to be valid 
for the strongly correlated electron systems. In the PAM, 
z^ 1 > z^ 1 is obtained from p T (0) < pj.(0) using eq. 
(11). The reason why the relationship Pj(0) < pj.(0) 
is derived can be explained as follows. From Eq. (3), 
At(0) = p 00 (—Ba — E CT (_B, u> = 0) + S (T (0,0)) is obtained 
within the DMFT at T = 0, where p 00 (w) is the DOS of 
/-electrons for U = and B = 0. In the present case, 
e/ + S(0) > p is satisfied (n = + n c < 2). We ap- 
proximate = Efe % + A* - e J - £«r (0, 0) - ^) 
considering that the w-dependence of DOS of conduction 
electrons is weak. p 00 (oj) increases (decreases) monoton- 
ically with ui for uj ~ in the case of n < 2 (n > 2). 
Thus, the relationship p-f (0) < pj. (0) is derived for B > 
at sufficiently low temperatures. Therefore the relation- 
ships z~ x > z" 1 and |ImS T (0)| > |ImSi(0)| are derived 
from eq. (11). 




Fig. 4. z T 1 /z l 1 (solid line), p x (0)/p T (0) (dotted line), and (1 - 

rt| )/(l — n^) (dashed and dotted line) against B in the 3D PAM 
for T = 0.01. 

In the following, we calculate spin-dependent mass en- 
hancements using the FLEX approximation and con- 
firm the validity of the analysis in the DMFT. In Fig. 
4, z-jT 1 1 ' z^ 1 (solid line), pj.(0)/pf(0) (dotted line) are 
plotted as a function of B in the 3D PAM for T = 
0.01. We confirm that z^ 1 / z^ 1 is almost proportional to 
pj,(0)/pf (0) in accordance with above relationship (11). 
In Fig. 4, (1 — n|)/(l — nj) (dashed and dotted line) is 
also shown against B, which represents z^ 1 / z^ 1 in the 
Gutzwiller approximation for U = oo. 15 As shown in Fig. 
4, (1 — nl)/(l — nf) is close to zZ 1 / z^ 1 obtained by the 
FLEX approximation in the PAM. Thus, z^ 1 / z^ 1 ob- 
tained by the FLEX and the Gutzwiller approximations 
satisfy the relationship (11). 



Finally, we discuss discrepancies between our obtained 
results and experiments. First one is that the obtained 
magnitude of z^ 1 is much smaller than that of the ex- 
perimental values with z^ 1 ^27 for majority spin and 
z^ 1 ~ 90 for minority spin, respectively, at H = 15T. 9 
Second one is that in our calculation the magnitude of 
z^ 1 increases and that of z^ 1 decreases with magnetic 
field, while in the experiment both z^ 1 and z^ 1 decrease 
with magnetic field. Recently, we found that this discrep- 
ancy may originate from the absence of vertex corrections 
for the susceptibility in the present calculations. 

We have shown that both Xzz and x^ decrease with 

magnetic field when the vertex corrections for suscepti- 
bility are included in the Hubbard model. 26 This result 
leads to the decrease in the magnitude of the mass en- 
hancement factors of both spins. 

In summary, we have obtained the mass enhancement 
factor (z^ 1 ) in the presence of the magnetic field for the 
PAM using the FLEX approximation. The resulting z^ 1 
is almost 10 and is much larger than that in the Hub- 
bard model. The ratio of the mass enhancement factor 
-i 

~~y increases with the magnetic field. Based on the per- 
turbation theory up to infinite order, we have derived a 
general relationships, i.e., ^Lt ~ by applying the 

DMFT. In the case of n < 2 (e f + E(0) > n), which cor- 
responds to Ce compounds, the relationship z^ 1 > z^ 1 
is obtained, since the DOS of the up-spin electrons at 
the Fermi energy is smaller than that of down-spin elec- 
trons. This result is consistent with the experimental 
data of de Haas- van Alphcn 9 in CeCoIns. We predict 
that z^ 1 < z^ 1 in Yb compounds, which correspond to 
n> 2 (e^ + S(0) < p). 
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